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The Stellar Activity - Rotation Relationship 
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Using a new catalog of 824 solar and late-type stars with X-ray luminosities and rotation periods we have studied the 
relationship between rotation and stellar activity. From an unbiased subset of this sample the power law slope of the 
unsaturated regime, Lx/Lboi oc Ro 13 , is fit as [3 = —2.70 ± 0.13. This is inconsistent with the canonical /3 = —2 
slope to a confidence of 5a and argues for an interface-type dynamo. Super-saturation is observed for the fastest rotators 
in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation 
radius and the excess polar updraft, the latter theory being supported by other observations. We also present a new X-ray 
population synthesis model of the mature stellar component of our Galaxy and use it to reproduce deep observations of a 
high Galactic latitude field. The model, XStar, can be used to test models of stellar spin-down and dynamo decay, as well 
as for estimating stellar X-ray contamination rates for non-stellar studies. 



1 Introduction 

Stars across the Hertzsprung-Russell diagram are known to 
emit X-rays with only a few exceptions. Solar- and late-type 
stars are thought to generate their X-rays from a magnet- 
ically confined plasma known as a corona (Vaiana et al. 
1981), which is believed to be driven by the stellar mag- 
netic dynamo, which itself is thought to be driven by dif- 
ferential rotation in the stellar interior (e.g., Parker 1955), a 
phenomenon that has been confirmed in the Sun through he- 
lioseismology (Duvall et al. 1984). The observed decrease 
in X-ray emission from the pre-main sequence to the age of 
the Galactic field population can therefore be attributed to 
the rotational spin-down of the star, which is driven by mass 
loss through a magnetized stellar wind (Skumanich 1972). 

The relationship between stellar rotation and tracers of 
magnetic activity is a vital probe of the stellar dynamo. A 
relationship between rotation and activity was first quanti- 
fied by Pallavicini et al. (1981), who found that X-ray lu- 
minosity scaled as Lx oc (vsmi) 19 , providing the first 
evidence for the dynamo-induced nature of stellar coronal 
activity. For very fast rotators the relationship was found 
to break down with X-ray luminosity reaching a saturation 
level of Lx/L bo i ~ 1CT 3 (Micela et al. 1985), independent 
of spectral type. This saturation level is reached at a rotation 
period that increases toward later spectral types (Pizzolato 
et al. 2003), but it is unclear what causes this. Despite much 
work there is yet to be a satisfactory dynamo theory that 
can explain both the solar dynamo and that of rapidly rotat- 
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ing stars and the continued lack of a sufficiently large and 
unbiased sample has no doubt contributed to this. 

We have produced a new catalog of stars with stel- 
lar rotation periods and X-ray luminosities, the details of 
which are presented in Wright et al. (2011). The cata- 
log includes 824 solar- and late-type stars, including 445 
field stars and 379 stars in nearby open clusters (ages 40- 
700 Myrs). The sample was homogenized by recalculat- 
ing all X-ray luminosities and converting them all onto the 
ROSAT 0.1 — 2.4 keV band. To minimize biases we re- 
moved all sources known to be X-ray variable, those that 
exhibit signs of accretion, or those in close binary systems. 
The sample is approximately equally distributed across the 
color range V - K s = 1.5 - 5.0 (G2 to M4) with —30 stars 
per subtype, dropping to —10 stars per subtype outside of 
this, from F7 to M6. 

Here we highlight some of the results derived from this 
sample, focussing attention on the regimes of low-activity 
dynamo decay and high-activity supersaturation. Details of 
these findings and further results and discussions from this 
sample can be found in Wright et al. (2011). Finally we 
present a new population synthesis model that can be used 
to address similar problems in stellar activity. 

2 Dynamo Efficiency in the Linear Regime 

The compiled catalog has been used to study the relation- 
ship between the level of X-ray activity, represented by the 
ratio of X-ray to bolometric luminosity, R\ = Lx / Lboi 
and the rotation period, parameterized by the spectral-type 
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Fig. 1 Rx = Lx /Lboi plotted against P rot (left) and the Rossby number, Ro = P rot jr (right), for all stars in our 
sample. Known binary stars are shown as plus symbols, and the Sun with a solar symbol. The best-fitting saturated and 
non-saturated activity-rotation relations are shown as a dashed red line in the right-hand panel. 



independent Rossby number, Ro = P rot /r, the ratio of the 
rotation period to the convective turnover time (Noyes et al. 
1984). Figure Q] shows Rx as a function of P rot and Ro, 
illustrating that the use of the latter parameter greatly re- 
duces the scatter in the rotation-activity diagram, and high- 
lights the two main regimes of coronal activity: a linear 
regime where activity increases with decreasing rotation 
period, and a saturated regime where the X-ray luminos- 
ity ratio is constant with log Rx = —3.13 ±0.08. The 
transition between these two regimes is found to occur at 
Ro = 0.13 ± 0.02 from a two-part power-law fit. 

The sample used here suffers from a number of biases, 
most importantly an X-ray luminosity bias due to the selec- 
tion only of sources with measured X-ray fluxes and pho- 
tometric rotation periods. To overcome this we used an X- 
ray unbiased subset of our sample, the 36 Mt. Wilson stars 
with measured rotation periods (Donahue et al. 1996), all of 
which are detected in X-rays. While some biases may still 
exist due to the ability to measure rotation periods, Don- 
ahue et al. (1996) conclude that any such biases are unlikely 
to affect the rotation period distribution. 
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Fig. 2 Rx = Lx/Lboi versus Rossby number for an 
unbiased sample of 36 stars with unsaturated X-ray emis- 
sion. The log-log OLS bisector fit, (3 — —2.70, is shown 
as a dashed line alongside a fit with the canonical slope of 
f3 = -2.0. 



A single-part power-law fit in the linear regime (where 
Rx oc Ro") to this sample is shown in Figure |2] with a 
fit of /3 = —2.70 ± 0.13. This is a steeper slope than fit 
using the full sample and notably steeper than the canoni- 
cal value of j3 — — 2 (Pallavicini et al. 1981), though their 
use of projected rotational velocities instead of rotation pe- 
riods represents a different relationship than that fitted here. 
Our slope is inconsistent with the canonical value to a con- 
fidence of 5cr, which argues against a distributed dynamo 
operating throughout the convection zone, the efficiency of 
which scales as Ro~ 2 (Noyes et al. 1984), and instead ar- 
gues for an interface dynamo (e.g., Parker 1993), which has 
a more complex dependency where (3^—2. 



3 Supersaturation 

At very high rotational velocities the fractional X-ray lu- 
minosity is observed to decrease below the saturation level 
(Randich et al. 1996), an effect dubbed supersaturation, 
the cause of which is still heavily debated. Figure [3] shows 
Lx/Lboi for all stars with saturated X-ray emission as a 
function of various quantities suggested as parameters of 
supersaturation. The first two panels show the effect of ro- 
tation period and Rossby number on the X-ray luminosity 
ratio, neither of which show evidence for supersaturation 
(i.e. a decline in Lx/L^oi) for either the fastest rotators or 
the sources with the smallest Rossby numbers. 
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Fig. 3 L \ I Lb i for stars with saturated X-ray emission as a function of the rotation period (upper left), the Rossby 
number (upper right), the Keplerian co-rotation radius (lower left), and the excess polar updraft (lower right). Stars are 
colored as per their spectral type: F-type (blue), G-type (green), K-type (black), and M-type (red). Also shown are 25-star 
average luminosity ratios, with standard errors. 



Two theories suggested to explain the supersaturation 
effect are coronal stripping and convective updrafts. In the 
theory of coronal stripping the star is spinning so fast that 
centrifugal forces strip away the outer layers of the corona, 
reducing the X-ray emitting volume (Jardine & Unruh 1999, 
James et al. 2000). This theory is parameterized by the coro- 
tation radius, Rxepier, the radius at which centrifugal forces 
balance gravitational forces in the corona. In Figure [3] it is 
clear that this parameter exhibits a correlation with the su- 
persaturation effect with the fastest rotators showing a drop 
to log R x = -3.41 ± 0.26. 

In the theory of convective updrafts, nonuniform heat- 
ing of the convective envelope (in accordance with the von 
Zeipel, 1924, theorem for rapidly rotating stars) results in a 
poleward migration of active regions and a reduction in the 
filling factor of active regions on the stellar surface (St§pieh 
et al. 2001). This can be quantified in terms of the excess 
polar updraft at the bottom of the convective envelope and 
is also shown in Figure[5J where the parameter also displays 



a strong correlation with the supersaturation effect with the 
fastest rotators dropping to log Rx = —3.45 ± 0.16. 

Based on this data it is difficult to distinguish between 
these two theories, partly because they have similar depen- 
dencies in terms of stellar parameters. An alternative dis- 
tinguishing method is the underlying coronal structure that 
the two theories require to be effective. Coronal stripping 
will occur if the loops are large, low density, and unstable to 
the centrifugal force. In Figure[3]a decrease in X-ray lumi- 
nosity begins at Rnepier/R* ~ 3, suggesting that coronal 
loops are as large as two stellar radii in height. Alternatively 
coronal loops may be compact, high density and influenced 
by the dynamics of surface flows where convective updrafts 
may influence the surface activity level. The difference be- 
tween the two supersaturation theories is therefore the size 
and density of the coronal plasma. High-resolution X-ray 
spectroscopy can be used to infer coronal plasma densities 
(and therefore scale heights), but such data for supersatu- 
rated stars is rare. The only such observation to date is of 
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Fig. 4 Cumulative distribution of X-ray fluxes for sources 
in the COSMOS field compared to two predictions from the 
XStar model for j3 = —2.0 and j3 = —2.7. The observations 
are complete to fx ~ 2 x 1CP 15 erg s _1 cm -2 . 



VW Cephei, for which coronal loop heights of 0.06-0.2 i?* 
have been measured (Huenemoerder et al. 2006). The rota- 
tional period of VW Cephei implies Rxepier / 'R* ~ 1.75, 
significantly higher than the observed coronal loop heights, 
and suggesting that coronal stripping is unlikely to be re- 
sponsible for supersaturation. 

4 Modeling the Decay of Stellar Dynamos 

Studies of the rotation - activity relation in the manner per- 
formed above are limited both by the sample size of stars 
with measured X-ray luminosities and rotation periods and 
also by biases inherent in compiling such samples. A dif- 
ferent approach that allows one to overcome such problems 
is to model the X-ray emission of a population of stars and 
compare the number and properties of such stars with obser- 
vations. We have compiled such a model, XStar, combin- 
ing a Galactic population synthesis model with models of 
rotational spin-down (e.g., Barnes 2003) and the activity- 
rotation framework outlined here. The results of this model 
can then be compared to X-ray observations that sample dif- 
ferent sight-lines in our Galaxy such as deep targeted obser- 
vations or wide-field, shallow surveys. 

An example of the model is shown in Figure |4] where 
two realizations of the model using (5 = —2.0 and j3 = 
—2.7 are compared to the cumulative flux distribution of 
X-ray sources in the high Galactic latitude Chandra COS- 
MOS field (Wright, Drake & Civano 2010). The models re- 
produce the observations very well down to the estimated 
observational completeness limit of fx ~ 2 x 10 -15 erg 
s _1 cm~ 2 , with similar forms and the brightest sources re- 
produced well. The j3 = —2.7 model provides the better fit, 
although various other parameters within the model must be 
fully tested before firm conclusions can be drawn. Nonethe- 
less the code does will in reproducing the properties (dis- 



tances, spectral types, X-ray fluxes) of the observations, and 
will therefore prove useful for studying the properties of 
mature, X-ray emitting stellar populations as well as for es- 
timating stellar contamination rates for studies of star form- 
ing regions and extra-Galactic populations (e.g., Brassing- 
tonetal. 2012). 

5 Summary 

A new catalog of stars with measured X-ray luminosities 
and rotation periods is used to study the rotation - activ- 
ity relation. The power-law slope of the unsaturated regime 
is fit as (3 = — 2.7 ± 0.13, inconsistent with the canoni- 
cal j3 = — 2 value with a confidence of 5<r and arguing for 
an interface-type dynamo. Coronal supersaturation in ultra- 
fast rotators is shown to correlate well with both the coro- 
tation radius and the excess polar updraft as the parameters 
of the coronal stripping and convective updrafts theories, 
though other observations support the latter theory. Finally 
a new population synthesis model of X-ray emitting stars in 
a Galactic sightline is introduced and tested successfully on 
deep observations of the high Galactic latitude COSMOS 
field. A better fit is found for /3 = -2.7 than for ft = -2, 
confirming results with a fixed sample. 
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